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Abstract 

Wc consider a directed graph on the 2-dimcnsional integer lattice, placing a directed 
edge from vertex (ii,^) to (31,32), whenever i\ < j±, ii < 32, with probability p, 
independently for each such pair of vertices. Let L n ^ m denote the maximum length of 
all paths contained in an n x m rectangle. We show that there is a positive exponent a, 
such that, if m/n a — > 1, as n — > oo, then a properly centered/rescaled version of L n>m 
converges weakly to the Tracy-Widom distribution. A generalization to graphs with 
non-constant probabilities is also discussed. 

Keywords and phrases. Random graph, last passage percolation, strong approxima- 
tion, Tracy-Widom distribution. 
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1 Introduction 



Random directed graphs form a class of stochastic models with applications in computer 
science [11], biology [5, 15, 16] and physics [10]. Perhaps the simplest of all such graphs 
is a directed version of the standard Erdos-Renyi random graph [2] on n vertices, defined 
as follows: For each pair {i,j} of distinct positive integers less than n, toss a coin with 
probability of head equal to p, < p < 1, independently from pair to pair; if head shows up 
then introduce an edge directed from mm(i,j) to m&x(i,j). There is a natural extension of 
this graph to the whole of Z studied in detail in [7]. In particular, if we define the asymptotic 
growth rate C = C(p), as the a.s. limit of the maximum length of all paths between 1 and 
n divided by n, [7] provides sharp bounds on C(p) for all values of p £ (0, 1). 

A natural generalization arises when we replace the total order of the vertex set by a 
partial order, usually implied by the structure of the vertex set. In such a model, coins are 
tossed only for pairs of vertices which are comparable in this partial order. The canonical 
case is to consider, as a vertex set, the 2-dimensional integer lattice Z x Z, equipped with 
the standard component- wise partial order: (11,12) H J2) if the two pairs are distinct 
and it < 12, ji < J2- Such a graph was considered in [6]. In that paper, it was shown [6, 
Thm. 5] that if L n t m denotes the maximum length of all paths of the graph, restricted to 
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{0, . . . , n} x {1, . . . , m}, then there is a positive k (depending on p and the fixed integer m), 
such that 

' L ^~ Cn \ t >o) (Z t , m , t>0), (1) 



K\ n 



where 2T. m is the stochastic process defined in terms of m independent standard Brownian 
motions, B^\ . . . , B^ m \ via the formula 



Z t , m := sup ^[Bf'-Bji], t>0. 

0=to<tl—<tm-l<*ra=t -_j 

One can speak of Z as a Brownian directed percolation model, the terminology stemming 
from the picture of a "weighted graph" on R x {1, . . . , m} where the weight of a segment 
[s, t] x {j} equals the change B^ — Bg of a Brownian motion. If a path from (0, 0) to (i, m) 
is defined as a union UjLi^i-i'^j] x {.?'} °f suc h segments, then Z represents the maximum 
weight of all such paths. 

Baryshnikov [3], answering an open question by Glynn and Whitt [8], showed that 

7 @ \ 

where A m is the largest eigenvalue of a GUE matrix of dimension m. Since Z. )Tn is 1 /2-self- 
similar, we see that 

7 w K\ 

Now, fluctuations of A m around the centering sequence have been quantified by Tracy 
and Widom [17] who showed the existence of a limiting law, denoted by Ftw : 

m 1/6 (A m - 2^m) F TW - 

m— >oo 

A natural question then, raised in [6], is whether one can obtain F^w 

as a weak limit of L n m 

when n and m tend to infinity simultaneously. Our paper is concerned with resolving this 
question. To see what scaling we can expect, rewrite the last display, for arbitrary t > 0, as 

mV 6 (% - 2^) Ftw . 
A statement of the form X(t,m) — - - > X, where the distribution of X(t,m) does not 

m— yoo 

depend on the choice of t > 0, implies the statement X(t,m(t)) — > X, for any function 

t— >oo 

mit) such that m(t) > oo. Hence, upon setting m = [t a ], we have 

t— >oo 

t^(^P--2VtA^F TW . (2) 

Therefore, it is reasonable to assume that an analogous limit theorem holds for a centered 
scaled version of the largest length L n r n ai, namely that 

( Ln ^ ] ~ Cin - 2yfc) F TW , (3) 
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where c\, C2 are appropriate constants. Since we are talking about interchange of limits here, 
it is also reasonable to assume that (3) holds provided that a is small enough. 

A different model which, a priori, seems to bear little resemblance to ours, is the directed 
last passage percolation model on Z d (the case d = 2 being of interest here) . We are given a 
collection of i.i.d. random variables indexed by elements of Zl. A path from the origin to 
the point n 6 Zi is a sequence of elements of Zi, starting from the origin and ending at n, 
such that the difference of successive members of the sequence is equal to the unit vector 
in the i-th direction, for some 1 < i < d. The weight of a path is the sum of the random 
variables associated with its members. Specializing to d = 2, let L n ^ m be the largest weight 
of all paths from (0, 0) to (n, m). Assuming that the random variables have a finite moment 
of order larger than 2, Bodineau and Martin [4] showed that (3) holds for all sufficiently 
small positive a (the threshold depending on the order of the finite moment). Independently, 
Baik and Suidan [ 1 ] obtained the same result for random variables with a finite 4-th moment 
and for a < 3/14. In both papers, partial sums of i.i.d. were approximated with Brownian 
motions, in the first case using the Komlos-Major-Tusnady (KMT) construction, while in 
the second using Skorokhod embedding. 

To show that (3) holds for our model, we adopt the technique introduced in [G], which is 
existence of skeleton points on each line Z x {j}. Skeleton points are, by definition, random 
points which are connected with all the other points on the same line. In [G] Denisov, Foss 
and Konstantopoulos used this fact, together with the fact that, for finite m, one can pick 
skeleton points common to all m lines, in order to prove (1). However, when m tends to 
infinity simultaneously with n, it is not possible to pick skeleton points common to all lines. 
Modifying the definition of skeleton points enables us to give a new proof of (1), as well as to 
prove (3). To achieve the latter, we borrow the idea of KMT coupling from [ :]. However, we 
need to do some work in order to express the random variable L n ^ m in a way that resembles 
a maximum of partial sums. 

Although we focus on the case where the edge probability p is constant, it is possible to 
consider a more general case, where the probability that a vertex (11,12) £ Z x Z connects 
to a vertex (ji, j'2) depends on the distances \ji — i\\ and \j2 — *2 1 of the two vertices. The 
ideas of the proof are the same. This generalization is discussed in the last section of the 
article. 



2 The one-dimensional directed random graph 

We summarize below some properties of the directed Erdos-Renyi graph on Z with con- 
nectivity probability p taken from [7]. For i < j, let be the maximum length of 
all paths with start and end points in the interval Then, for i < j < k, we have 
L[i, k] < L[i,j] + L[j, k] + 1. Since the distribution of the random graph is invariant under 
translations, and is also ergodic (the natural invariant cr-field is trivial), it follows, from 
Kingman's subadditive ergodic theorem, that there is a deterministic constant C = C(p) 
such that 

lim L[l,n]/n = C, a.s. (4) 

n— >oo 

In fact, C = inf n >i EL[1, n]/n. The function C(p) is not known explicitly; only bounds are 
known [7, Thm. 10.1]. For example, 0.5679 < C(l/2) < 0.5961. We also know that there 
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exists, almost surely, a random integer sequence {T r ,r E Z} with the property that for all 
r, all i < T r , and all j > T r , there is a path from i to T r and a path from T r to j. The 
existence of such points, referred to as skeleton points, is not hard to establish ([6]). Since 
the directed Erdos-Renyi graph is invariant under translations, so is the sequence of skeleton 
points, i.e., {T r ,r E Z} has the same law as {n + T r ,r E Z}, for all n E Z. Moreover, it 
turns out that the sequence forms a stationary renewal process. If we enumerate the skeleton 
points according to • • • < T_i < Tq < < Ti < ■ ■ ■ , we have that {r r+ i — T r , r E Z} are 
independent random variables, whereas {r r+ i — T r ,r ^ 0} are i.i.d. Stationarity implies 
that the law of the omitted difference T\ — Tq has a density which is proportional to the 
tail of the distribution of T2 — Ti. In [6] it is shown that the distance T 2 — T\ between 
two successive skeleton points has a finite 2nd moment. One can follow the same steps of 
the proof, to show that in our case, with constant probability p, this random variable has 
moments of all orders. Moreover, one can show that for some a > (the maximal such a 
depends on p) it holds that Ee a ^- 2 ~ Tl ^ < 00. 

The rate Ao of the sequence of skeleton points can be expressed as an infinite product: 

1 00 

For example, for p = 1/2, Ao ~ 1/12. 

A central limit theorem for L[l,n] is also available [6, Thm. 2]. If we let 

(rg-var^^r^-C^-ri)), (6) 

then 



where N(0, 1) is a standard normal random variable. Note that <Tq 7^ var(L[ri, I^]). Unfor- 



tunately, we have no estimates for Oq, but, interestingly, there is a technique for estimating 
it, based on perfect simulation. This was briefly explained in [7] in connection with an 
infinite-dimensional Markov chain which carries most of the information about the law of 
the directed Erdos-Renyi random graph. 

In addition, it is shown in [7] that C can also be expressed as 

EL[T U T 2 ] 

c -E(T 2 -r 1 y (8) 

In fact, if \y r ,r E Z} is a random sequence of integers, defined on the same probability 
space as the one supporting the random graph, such that {T Ur , r E Z} is a stationary point 
process then 

_ EL[T Vr ,T Ur+1 ] 



£(rv +1 -r 



V r 



The most important property of the skeleton points is that if 7 is a skeleton point, and 
if i < 7 < j, then a path with length L[i,j] (a maximum length path) must necessarily 
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contain 7. This crucial property will be used several times below, especially since, for every 
i < j, the following equality holds 

L\Ti,Tj] = L[T i ,T i+1 ] + L[T i+1 ,T i+2 ] + --- + L[T j ^ l ,T j }. 

Furthermore, the restriction of the graph on the interval between two successive skeleton 
points is independent of the restriction on the complement of the interval; hence the sum- 
mands in the right-hand side of the last display are independent random variables. 



3 Statement of the main result 



It is clear from (7) that the constants c\,c 2 in (3) should be as follows: c\ = C, c 2 = V Act 2 . 
Now we can formulate the main result. 

Theorem 1. Let C , Ao, ct 2 be the quantities associated with the directed random graph on 
Z with connectivity probability p, defined by (4) (equivalently, (8) ), (5), (6), respectively. 
Consider the directed random graph on Z x Z and let L„ )m be the maximum length of all 
paths between two vertices in [0, n] x [l,m]. Then, for all < a < 3/14, 

(9) 

where F?w i> s the Tracy-Widom distribution. 



To prove this theorem, we will first define the notion of skeleton points for the graph on 
Z x Z and then prove pathwise upper and lower bounds for L nm which depend on paths 
going through these skeleton points. This will be done in Section 4. In Section 5.1 we show 
that the difference between these bounds is of the order o(n b ), where b = (1/2) — (a/6) 
is the net exponent in the denominator of (9). We will then (Section 5.2) introduce a 
quantity S n ,m which resembles a last passage percolation problem and show that it differs 
from L n ^ m by a quantity which is of the order o(n b ), when m = [n a \. The problem will then 
be translated to a last passage percolation problem (with the exception of random indices) . 
This will finally, in Section 5.3 be compared to the Brownian directed percolation problem 
by means of strong coupling. 



4 Skeleton points and pathwise bounds 



Our model is a directed random graph G with vertices Z x Z. For each pair of vertices i, 
j, such that i -i j, toss an independent coin with probability of heads equal to p; if a head 
shows up introduce an edge directed from i to j. 

A path of length I in the graph is a sequence (*0j *i, • • • ; it) of vertices io -< ii -< ■ ■ ■ -< %i 
such that there is an edge between any consecutive vertices. 
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We denote by G n>m the restriction of G on the set of vertices {0,1, ... ,n} x {1, . . . , m}. 
The random variable of interest is 

L n ,m '■= the maximum length of all paths in G n ^ m . 

We refer to the set Z x {j} as "line j" or "j-th line", and note that the restriction of G 
onto Z x {j} is a directed Erdos-Renyi random graph. We denote this restriction by G^\ 
Typically, a superscript (j) will refer to a quantity associated with this restriction. For 
example, for a < b, 

L^[a,b] := the maximum length of all paths in G^' with vertices between (a,j) and (b,j) 
and we agree that L^[a, b] = if a > b. 

Clearly, the {G^\j e Z} are i.i.d. random graphs, identical in distribution to the directed 
Erdos-Renyi random graph. Therefore, for each j £ Z, 

lim L^\\, n] In = C, a.s. 

n— >oo 

To establish upper and lower bounds for L njTn , we need to slightly change the definition 
of a skeleton point in G. 

Definition 1 (Skeleton points in G). A vertex of the directed random graph G is called 
skeleton point if it is a skeleton point for G^> (for any i' < i < i" , there is a path from (i' ,j) 
to (i,j) and a path from (i,j) to (i",j)) and if there is an edge from (i,j) to + 1). 



Therefore, the skeleton points on line j are obtained from the skeleton point sequence of 
the directed Erdos-Renyi random graph G^' by independent thinning with probability p. 
When we refer to skeleton points on line j, we shall be speaking of this thinned sequence. 
The elements of this sequence are denoted by 

• • • < r°] < r£' } < o < if } < if > < • • • 



and have rate 

1 



v^=p\\{i-{i-p) k f. 



£ J \ L 2 ~ L i ) k=l 
The associated counting process of skeleton points on line j is defined by 

$W(f)-#(s) = J]l(s<rW<t), s,teR, s<t, 
together with the agreement that 

= 0. 

Note that we insist on having the parameter t in <fr(i\t) as an element of R (and not just 
Z). We also let 

x(o)( f \ — r0') 

y0)(t) : = r^ (i)(t)+1 , 
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be the skeleton points on line j straddling t: 

X&(t)<t<Y®(t). (10) 

Next we prove upper and lower bounds for L n;W ,. The set of dissections of the interval 
[0, n] Clinm non-overlapping, possibly empty intervals is denoted by 

T n ,m '■= {t = (to, h,..., t m ) £ M m+1 : = to < ti < ■ ■ ■ < t m ^i < t m = n}. 

Lemma 1. (Upper bound) Define 

m 

I n , m := sup ^L^lX^it^), Y^(tj)]+m. (11) 

Then L nTn — -^n,m* 

Proof. Let ir be a path in G n ,m- Consider the lines visited by ir, denoting their indices by 

1 < v\ < U2 < • • • < vj < m. Let (oj, Vj) and (bj, Vj) be the first and the last vertex of line 
Vj in the path ir. Then the length of ir satisfies 

J 

Since successive vertices in the path should be increasing in the order -<!, we have < aj, 

2 < j < J. Hence, with bo ■= 0, 

J J 
M <J2 L{Vj) i b 3-ii b j} + J~1< ^LMixWQij-!), Y^(bj)} + J - 1, 
i=i i=i 

where we used (10). Since J < m, we can extend O = 6o<6i<---<foj<ntoa dissection 
of [0, n] into m non-overlapping intervals, showing that the right-hand side of the last display 
is bounded above by L> n ,m- Taking the maximum over all 7r in G n ^ m , we obtain L n ^ m < L n ^ m , 
as required. □ 

Note that the existence and properties of skeleton points were not used in the proof of 
the upper bound, other than to ensure that the upper bound is a.s. finite. 

Lemma 2. (Lower bound) Define 

Ag> := max (T& - if), 

0<i<*W(n) 



and 



L njm := sup gLWfKW^), ^(t^-^Att). 



ThcTl L n ^ m > L.n,m' 
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Proof. We will show that, for all t = (to, • • • ,t n ) £ 7~ n ,m, there is a path tt in G nyTn with 
length \tt\ satisfying 

m m 

E^'^fe-l). XU) ^)] < Kl + E A n- (12) 
3=1 3=1 

Fix t G 7n m an d use the notation 

/ 3 = [^fe-i). ^ a) fe)] = [«3^3] ; J = l,...,m. 

Note that aj > 6j if there is one or no skeleton points on the segment (tj—\,tj] x {j} and 
then L^\lj) = 0. 

Given two skeleton points (x,i), (y,j) we say that there is a staircase path from (x,i) to 
(y,j) if there is a sequence of skeleton points 

(x,i) = (x ,i), (xi,i + l),..., (xj-i,j) = (y,j), 

such that x = xq < x\ < ■ ■ ■ < Xj_i = y. See Figure 1. Clearly then, there is a path from 
(x,i) to (y,j) which jumps upwards by one step each time it meets a new skeleton point 
from the sequence. We denote this by 

(x,i) -i, (y,j). 



(j) 



-X X-X- 



(O 





1 


\ 

X > 


\ 

1 

< x 1 — 



Figure 1: A staircase path from (x,i) to (y,j) jumps upwards at skeleton points (denoted by 
x) but may skip several of them before deciding to make a jump 



Among all the staircase paths from (x,i) to (y,j), we will consider the best one, defined 
by two properties: 

• Property 1: A best path from (x,i) to (y,j) jumps from line k to line k + 1, k = 
i, i + 1, . . . ,j — 1, at the first next skeleton point on line k, i.e. at the points xo and 
x k _ i+l = Y( k+1 Xx k -i), k = i,...,j -2. 

• Property 2: Every horizontal segment of a best path is a path of maximal length. 
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If all the intervals I\, . . . , I m are empty, the left-hand side of (12) is zero and the inequality 
is trivially satisfied for any path tt. 

Otherwise, for a fixed t £ G n ^ m we will construct a path tt in G n , m for which (12) holds. 
Define a subsequence v\ < zv 2 < • • • of 1, . . . , m, inductively, as follows: 

i^i := inf{l <j<m: Ij ^ 0}, (13) 
v r := inf{j > i/ r _i : (6 I , T ._- L , tV-i) A r > 2. (14) 

See Figure 2 for an illustration. The procedure stops if one of the elements of the subsequence 
exceeds 772 or if the condition inside the infiniuni is not satisfied by a path in Gn,m- 

(9) - 

a g c 8 b g 

(8) 
(7) 



(5) 
(4) 
(3) 
(2) 
(1) 



x 


Xb 


) 


— x- 

l 










I 




c 5 a 5 b 5 






— X — X X— 




V 

\s \ ' 


^ a. 
X 










b 3 ' 




* 
















7\ 




-x — 






— x — 
















X 
























XX 





Figure 2: Illustration of the procedure defined by (13)-(14). There are four best staircase 
paths: the path from (62,2) to (63,3), the path from (63,3) to (64, 4), the path from (64, 4) 
to (65,5), and the path from (£>5,5) to (b$,8). Observe that Ij = (aj,bj), in the figure, are 
nonempty only for j = 2, 5, 7 and 8 (these are the highlighted intervals), but I7 is not visited 
by the constructed path. Moreover, Is is only partly visited and the path enters Is at a point 
cs between bs- 

Let J be the last index in the above defined sequence. Let tt\ be a path of maximum 
length from (a Vl , v\) to (b Vl , v\) and define, for r = 2, 3, . . . , J, a path u r as a best staircase 
path from {b Ur _ x , u r -i) to (b Ur ,u r ). Note that, for each r = 2,3, ...,J, J > 2, the end 
vertex of 7r r _i is the start vertex of TT r . Therefore we can concatenate the paths 7Ti, . . . , ttj 
to obtain a path tt. This path starts from (a Ul , u\) and ends at (b Vj , vj). 

Let 

tt^' := the restriction of path tt on line j 
and 1 7r^') I its length on line j. Also, for j > v\ denote 

(cj,j) := the first vertex on line j of path tt. 

Split the sum in the left-hand side of (12) along the elements of the subsequence {v\, . . . , uj}: 

m J+l i> T J+l 

j=l r=Xj=u r - 1 +l r=l 
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where we have conveniently set 

v Q : = 0, v J+1 := m, 

in order to take care of the first and last terms. By the defintion of u%, the intervals 
1%, I2, ■ ■ ■ , I vi -i are empty and 

G 1 =L^(I U1 ) = \tt^\. 
Assume now that 2 < r < J, and write 

v r v r — l 

G r = £ L®(I j ) = Y L^{I 3 ) + L^\l Ur ). 

j=v r -i+l j=v r -i+l 

Since ir Ur is the path of maximal length from (c Ur ,i/ r ) to its end- vertex (b Ur ,u r ) (Property 
2), if c Vr < a Ur then L^ Ur \l Ur ) < \ir^ Ur ^\. Define in this case l' v = 0. Otherwise, we can write 

L^\l Vr ) < L^[a Ur ,Y^\c Vr )] + L^[Y^\c Ur )X r ]. 

Then, again because of Property 2, L^\Y^ Vr \c Vr ), b Vr ) < \-k^\ and it is left to find a 
bound on the interval I' Ur = [a Ur , Y^ Ur \c Vr )]. Recall that by Property 1, depending whether 
j is a member of the sequence {v r ,r = 1,...,J} or not, Cj+i = bj or c J+ i = Y^\cj), 
respectively. Also, because of Ij C [tj-\,tj], we know that L^(Ij) < tj — ij-i- Hence, if 
v r — v r -\ > 1 it holds 

v r — 1 V r — 1 

j=fr-l+l j = ^ r _l + l 

I/ r V r 

<Y^(c Ur )-b Ur _ 1= &™ (<*)-<*)< E A n ] - 

j=v r -i+l j=V r -l+l 

Combining the above, we obtain 

v r — 1 

G r < Y A^ + |vr^)|. 

j = V r + l 

\ivj = m, then G j + \ = 0. Otherwise, we can extend the sequence {cj , j = v\ , v\ + 1, . . . , uj} 
defining iteratively c Vj+ \ := b V] and c J+ i := Y^\cj) until Cj > n for some j. Let if be the 
last index such that ck < As there was not possible to construct the best staircase path 
after the line z/j, K is at most m. Similarly as above, for Grj+i it holds 

m m 

G J+1 = Y ^V;)< Y (tj-tj-i)<n-b vj 

j = Vj + l j = vj + l 

= Y (Y {j \cj) - Cj ) < y 

j = Vj + l j = vj + l 



Finally, we obtain 

m m J m 

Yl U) (V < E A ^ } +Ei- K) i <E A ^ J) + m. 

j=l j=\ r =l j=l 

as required. □ 
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5 Further estimates in probability and Brownian directed 
percolation 



In the present section we prove Theorem 1 as a sequence of lemmas. 



5.1 Asymptotic coincidence of the two bounds 

Looking at (3), we can see that the correct scaling requires exponent 

2 6 

in the denominator and condition a < 3/7, which is equivalent to a < b. 

In the following two lemmas we will not specifically use the definiton of b and condition 
on a. Both lemmas hold for more general a, 6 > 0, < 6 — a < 1. 

Lemma 3. With b = (1/2) - (a/6) and a < 3/7, 

L n ,[n a ] ~ L. n ,[n a ] (p) 



> 0. 



Proof. Let t be such that the maximum in the right-hand side of (11) is achieved. Then 

mm m 

Ln,m - L n , m < m + Ag> + £ LW[JTW(t,-_i), - £ L^Wfe-i), X^fe) 

3=1 3=1 3=1 

m m s 

< m + A#> + I L^[xW (tj-!), rW(tj_i)] + ^ {i) (ii 

3=1 3=1 ^ 

m m 

< m + Y A £) + 2 V max L ^ [if } , r^J 

m 
3=1 



Hence 



n° L J ' l ' n° n° 



Since b > a and the random variables {Tv/j — I"^ 1 , i > 1} have a finite 1/(6 — a)-th moment, 
the converegence to for the second term above follows by Lemma A.l. □ 



5.2 Centering 

We introduce the quantity 

m f 

Sn,m ■= SUp J2 \ ^[X^itj-^X^itj)] ~ C[X^\t 3 ) - iWfe-i)] 

This should be "comparable" to L m ^ n — Cn when m = [n a ]. Indeed, we have: 
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Lemma 4. With b = (1/2) - (a/6), and a < 3/7, 

S n ,[n a ] ~ (L n ,[n a ] ~ Cn) (p) 



> 0. 



Proof. We begin by rewriting the numerator above as 

= snp \ Y,L {j \X^\t^ 1 ),X^\t 3 )\ + Cn - CY,[X {] \t 3 ) - X^Xt^l - L n)m . 

t67n,m L 7=1 j = 1 

Upon writing n = YljLiih ~ f° r an y * e 7n,mj we have 

m mm m 

j=i j=i j=i j=i 

Hence, on the one hand we have 

m m 

S n ,m ~ (L n , m ~ Cn) < sup LU) l X ® X0) - L ™>™ + 2C E A « ° 



*6T„. 



'3=1 



< sup ^LW[xW(t i _ 1 ),yW(t i )]-^,m + 2C^AW 

m 

<In,m-i„,m + 2C^AW. 

3=1 



On the other hand, 

m m 
S„, m " (L„,m - Cn) > sup £ [X & (tj^), X ® (tj)] - L n , m - 2C£ A#> 



teTn 



3=1 



3=1 



> sup ^i^iyW^o.iWfe)] - L n , m - 



3=1 



3=1 



> L„ j?n - L„, m - 2C £ ^ = -(Ln,m - L n) J - 2C ]T ■ 

3=1 3=1 

Therefore, 

m 

\S n>m - (L n>m -Cn)\< L n>m - L n , m + 2CJ2^ ] - 

3=1 

Thus, for m = [n a ], the result follows by applying Lemma A.l and Lemma 3. 
Define now variance a 2 as 

a 2 :=var(LW[r£ 1 ,ra-c'(r? ) -r£ 1 )) 



□ 
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and observe that a 2 = a 2 /p. We work with the quantity ^S n , m , which can be rewritten as 

-s n , m = sup Yl x< k ] > 

teTn ' m i=ifc=$W(t ; ,_ 1 )+i 

where 

Note that the random variables {xi"^}fc>i,j>i) indexed by both k and j, are independent 
and that {Xk^}k>2.j>i are identically distributed with zero mean and unit variance. The 
fact that the {x^}j>i do not have the same distribution will not affect the result, so we 
will not separately take care of it. 



5.3 Coupling with Brownian motion 

The term ^S n ^ m resembles a centered last passage percolation path weight [4], except that 
random indices are involved. Therefore, we start using the idea of strong coupling with 
Brownian motions, analogously to the proof in [4]. Let B^\B^ 2 \... be i.i.d. standard 
Brownian motions, and recall that 



m 

Z n , m := sup Y^-E^J. (15) 
Define the random walks , , . . . by 

^E*?, < = 0,1,2,..., (16) 



teT ri 



fe=i 



with Rq = 0. With this notation, we have 

1 m 

-S n , m = sup " R % tj J- ( 17 ) 

tfc In, in j = 1 

Taking into account (2) and Lemma 4, it is evident that to prove Theorem 1 it remains to 
show that 

a 1 £>n,[n a ] ~~ V^X-^n,[n a ] (p) 



11^ n— >co 



>0, 



or, using the scaling property of Brownian motion, that: 
Lemma 5. For all a < 3/14, 

° S n ,[n a ] ~~ ^Xn,[n a ] (p) 



> 0. 



To show that the random walks are close enough to the Brownian motion we use the 
following version of the Komlos-Major-Tusnady strong approximation result [12, Thm. 4]: 
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Theorem 2. For any 0<r<l, n£N and x £ [ci(log n) 1 / 1 ", C2(n log n) 1//2 ] , starting with a 
probability space supporting independent Brownian motions j = 1,2,..., we can jointly 
construct i.i.d. sequences = (Xi 1X2 ?■■■)> = 1 5 2, . . . , tuif/i £/je correct distributions 
and, moreover, such that, with pSp as in (16) above, 

P( max j R 0) - > x) < Cnexp{-ax r } /or a// j = 1, 2, . . . , 

l<i<n 

where the constants C,c\,C2 are depending on a, r and distributions ofxiP an d X2 • 

In addition to this, in order to take care of the random indices appearing in (17), we 
need a convergence rate result for the counting processes > 1} which is proven m 

the appendix. 



Proof of Lemma 5. From (17) and (15) we have 

m ( 1 

k"Xm " Z, n , m \ < sup £ \(B% - ~ (<. " B^J] 



teTn, 



■3=1 



oiit^ J I - rC?) I 1 I pO') _ R (i) I 



< 



+ \B [j) -B {j) \ + \B {j) - B {i) 



where 



3=1 3=1 

U<P := max - B<?\ V® := sup |^ - B#| 

0<i<n 0<s<n w 



Therefore, it is enough to show that, 

[n a ] [n a ] 

i-Vf/WJ^o and i-VyW^Uo. 

n b n-voo n b ' n-j-oo 

3=1 3=1 

For the first convergence we will take into account the coupling estimate as in Theorem 2. 
That is, we will throughout assume that the random walks and Brownian motions have been 
constructed jointly. The second convergence will be established without this estimate, i.e., 
we will show that it is true, regardless of the joint construction of the Brownian motions 
and the random walks. This is because the coupling we use is not detailed enough to give 
us information about the joint distribution of B^> and the counting process (the latter 
is not a function of the random walks used in the coupling). 



Proof of the first convergence. Let 8 > 0. We need to show that 

[«*"] 

P {J2 U n } > Snb ) ^° 

as n — > 00. 

3=1 
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Let e < b — a. Then 

[n a ] [n a ] 

P( V > 5n b ) < P( max < n e , V U® > 5n b ) + P( max U® > n £ ) . (18) 

' l<j'<[n a l — ' Kj<\n a ] 

j=l ~ J - 1 3 =l - J ~ l 

The first term from right-hand side is zero for large n. We are allowed, for n large enough, 
to estimate the second term using Theorem 2 for an arbitrary r £ (0, 1) and x = n £ as 

P( max Ui j) > n e ) < n a P( max \Bp - rP\ > n £ ) < n a Cnexp{-an er } -> 0, 

1<J<M /_ V l<i<n 1 

as n ^ oo. 

Proof of the second convergence. Let 1/4 < e < b — a. Replacing Un by Vn in (18), 
we see that the first term is again zero for large n and it remains to show that second term 
converge to 0, i.e., it is enough to show the convergence for its upper-bound 

n a P(V^ > n e ) -> 0, 

as n — > oo. Let 7 > 1 and 1/2 < q < 2e. Then we can write 

n a P(V^ > n e ) < n a P( sup |S«(s) - As| > 7 n 9 ) 

0<s<2n 

+ n a P( sup - > n e , sup |*«( 8 ) - Xs\ < 7 n"). 

0<s<n 1 ; 0<s<2n 

Corollary A.l implies the convergence of the first term above 

n a P( sup |$W(s) _ As] > 7n 9 ) < ra a exp{-a(2n) ?r } -> 0, 

0<s<2n 

as n — > 00, where 0<r<2 — 1/q. Set 99 = jn q /A. For the second term using the fact that, 
under our condition, As — 7n 9 < $( x )(s) < As + 7 n<? for s € [0, 2n], we have 



a P( sup \B% -B%\>n e , sup |$W( S )-A S | < 7 n«) 

0<s<n. 1 J 0<s<2n 

<n a P( max sup \b£L„ ,-B[!Lj>ti £ , sup (s) - Asl < >yn q ) 

< n a P( max sup IbQ , . ,, - B^, , J > n £ ) 

--yn1<t<'yn q 

< n a (n/ip + 1)P( sup I B t - B s \ > n e ) 

0<s<t<3jni 

< n a {n/ip + 1)P( sup \B S \ > n £ /2) 

0<s<3'yni 

< An a {n/ip + l)P{B- iinq > n £ /2) (a) 

< 4n a (n/<p + 1) exp{-n 2 7(24 7 n <? )} -»• 0, (b) 

as n — > 00. The inequality (a) is the consequence of P(sup 0<s<t B s > x) = 2P(B t > x) and 
the inequality (b) is an estimate for the tail of the normal distribution. □ 

Remark 1. The condition a < 3/14 is equivalent to b— a > 1/4 and, thus, it ensures existence 
of an e > 1/4 and later existance of 1/2 < q < 2e. Therefore, it is necessary for application 
of Lemma A.l, as well as for the convergence of the last estimate above. 
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6 Graph with non-constant edge probabilities 



In [6] the authors consider an one-dimensional model with connectivity probabilities de- 
pending on the distance between points, i.e. there is an edge between the vertices i and j 
with probability P\i-j\- To prove a central limit theorem for the maximal path length in 
that graph, two conditions are introduced: 



The conditions guarantee the existence of skeleton points and finite variance, defined as in 



We can extend the idea of the graph on Z x Z keeping, whenever possible, the notation 
from the Z x Z graph with constant probability p. 

Let (pij,i,j € No) be a sequence of probabilities that satisfies the following: 

• < pi, < 1; 

• < p ,l < 1; 

• Z)fcli fc r-1 (l -Pl,o)(l -P2,o) • • • (1 -Pk,o) < oo for some r > 2. 

Then the vertices (11,12) and (j'1,,7'2), (h^h) -< Oi, ^2)5 are connected with probability 

Pjl—il,j2—i2 • 

The conditions above are needed to ensure the existence of skeleton points, as in Definition 
1, and a finite r-th moment of the random variables X2, X3i • • • f° r some r > 2. 

From now on, let r denote the order of the highest finite moment, i.e. 



• < px < 1; 



• £r=iMi-pi)---(i-^)<oo. 



(6). 



r = 



sup{g > : kQ X (! -Pi,o)(l -P2,o) •••(!- Pk,o) < 00}. 



k=l 



Now we can state the analogue of Theorem 1 for this more general setting: 
Theorem 3. For all a < min(3/14, (r - 2)/(3r/7 + 1)), 
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Proof. We follow the lines of the proof of Theorem 1, emphasizing the points one should be 
careful about or which need to be modified. 

The construction of the upper and lower bounds for L n i n a] is the same as in Section 4. 
In Lemma 3 and Lemma 4, the term n a E[Ai]/n b converges to if a + 1/r < b, which leads 
to the constraint a < 6/7(1/2 — 1/r). 

Next, we rewrite, as in Lemma 5, 

m m 

\o- l S n , m - Z Xn>m \ < 2 + 2 V® 

3=1 3=1 

and prove convergence of each term separately. 

Proof of the first convergence. Instead of Theorem 2, we use the combination of results 
from Komlos-Major-Tusnady [12] and Major [14], as stated in Proposition 2 of [4], which 
allow us to couple Brownian motions j = 1, 2, . . . and random walks Ry\ j = 1, 2, . . . 

so that 

P( max \B& - i#| > x) < Cnx~ r , for all n G N, and all x G [n 1/r , n 1/2 }. (19) 

l<i<n 

Let e = 1/2. We want to establish convergence of the same terms as in (18). Applying, on 
the first term, the Markov inequality and properties of the coupling with Brownian motion 
(19) yield: 

P( i<m< ax C#) < n y\ I] U<p > Sn b ) < ^[11^1(11^ < n^)} 
i<j<[n j=1 n 

a ( rn 1 / 2 \ a ( rn 1 ? 2 \ 

5 h + L p(f, » ' =■ x)ix ) - h [ nllr + L Cnx ^ dx ) 

on" \ r — 1 r — 1 J \ r — 1 J on" 

as n — > oo. For the second term we use again the coupling properties (19): 

P{ max > n 1 / 2 ) < n a P( max - R a) | > n £ ) < n a Cnn- r ' 2 -> 0, 

l<i<["- a ] v l<i<n' * 1 1 

as 7i — \ oo because o + l — r/2<0 due to our new constraint on a. 

Proof of the second convergence. Here, the only change is the replacement of Lemma A. 2 
by part of Theorem 6.12.1 in [9], in our notation: 

If 1/2 < q < 1, then for all e > it holds that n 9r ~ 2 P(max < fc < n |r fc - Xk\ > en q ) ->■ as 
n — > oo. 

One can, in the same fashion as before, prove the following analogue of Corollary A.l, 
n qr ~ 2 P( sup |$(t) - At | > en q ) ->■ as n oo. 

0<t<n 

Likewise in Remark 1, condition a < 3/14 is necessary. Another constraint that occurs 
is a < qr — 2 and it can be shown that this is satisfied if a < (r — 2)/(3r/7 + 1). □ 
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A Appendix 



Lemma A.l. Let r > 1 and > 1} be a sequence of non-negative i.i.d. random 

variables such that E\Xi\ r < oo. Then, 



max Xi 

Ki<n 



as n — } oo. 



Proof. It suffices to prove the lemma for r = 1, since then the result for r > 1 follows easily 
by Jensen's inequality. Let M n = maxi<j< n Xi and S n = X\ +X2 + • • • +X n . Borel-Cantelli 
lemma, together with the assumption EX < 00, give X n /n — > 0, almost surely, as n — > 00. 
An easy computation shows that then also M n /n — > 0, almost surely, as n — > 00. 

On the other hand, we know that M n < S n for all n and that {S n /n, n > 1} is uniformly 
integrable. Therefore, M n /n is also unformly integrable and EM n /n — > as n — > 00 [9, 
Thm 4.5, Thm. 5.2]. □ 

The following lemma is an analogue of a result by Lanzinger [13, Prop. 2], with slight 
improvement in the probability bound. Moreover, we specialize the lemma to our case. 

Lemma A. 2. Let 1/2 < q < 1 and r < 2—1/q. Suppose that X, X%, X2, ■ ■ ■ are i.i.d. random 
variables such that EX = fi and Ee a \ x \ < 00 for some a > 0, and set S n = Y^7=i -^k- Then, 
for all e > 1, 

exp{an qr }P( max \Sk — kfi\ > en q ) — > 

as n — > 00. 



Proof. Without loss of generality we can assume that /i = 0. We first show that, for all 
e > 1, ex.p{an qr }P(S n > en q ) — > as n — > 00 . 

For fixed n G N, define X' k = X k l(X k < en q ) and set S' n = £jj =1 X' k . Note that EX' < 
and that for a', where a' < a and a'e T > a, it holds that EX 2 e a '^ x ^ r < 00. We can write 

P(S n > en q ) < P(S' n > en q ) + nP(X > en q ). (20) 

We first find a bound for the moment generating function for X' in a'(en q ) r ~ 1 : 

Ee a>{en*r-iX> < 1 + l a '\ £n ^(r~D {EX '2 + EX^W 1 *' l{X' > 0}) 

< i + l a ' 2 {en q ) 2( - r - l \EX 2 + EX 2 e a '\ x \ r ) 

< 1 + Cn 2 ^-^ < exp{C7n- 2g(1 - r) }, 

where for the first inequality we used e y < l + max{l, e y }y 2 /2 and for the last one 1+y < e y . 
Now, using Markov's inequality and the bound above for the first term in (20) yields 

P(S' n > en q ) = P ^(^r-'S' n > e «'K) r ) < e -«'(^r ^Ee a '( £nq y- lx 'y 

< exp{-c/ (en q ) r + Cn n - 2q{1 - r) }. 
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For the second term in (20), again using Markov's inequality, we have 

nP(X > en q ) = nP{e aXr > e a ^ r ) < ne- a( - £n ^ r E{e a \ x \ r ). 

Combining the two estimates finally establishes that 

exp{an qr }P(S n > en q ) < exp{-n qr (a' e r -a- Cn l - 2q+qr )} 

+ n^p{-n qr a{e r - l)}Ee a ^ r -> 

as n — > oo because of the choice of r, a' and e. 

By symmetry, the same convergence rate holds also for P(S n < —en q ). 

Thus, the statement of the theorem follows using Levy inequality (Theorem 3.7.2 in [9]) 
and the fact that for n large enough we can find e' such that 1 < e' < e — V2a 2 n l ^ 2 ~ q : 



P{ max |S fc | > en q ) < 2P{\S n \ > en q - V2na 2 ) = 2P(\S n \ > n q (e - si^n 1 ' 2 '* 1 )) 

l<k<n 

< 2P(\S n \ > e'n q ). 

□ 

Corollary A.l. Let X, X±, X2, ■ ■ ■ be positive, integer-valued i.i.d. random variables and 
suppose that the assumptions of Lemma A. 2 are satisfied. Then, for the counting process 
where $(f) = max{n : S n < t}, it holds that 

exp{an qr }P( sup - 1\ > en q ) ->■ for all e > 1. 

0<t<n 

Proof. For the counting process we have t — -X"$(t)+i < S$m < t. Therefore, we can write 
{ sup \/jL$(t) -t\> en q } = { sup (//*(<) - 1) > en q } U { inf (//*(*) -t)< -en q } 

0<t<n 0<t<n 0<t<n 

C { sup - S m ) > en q } U { inf - S m - X m+1 ) < -en 9 } 

0<t<n 0<t<n 

C { sup (//$(t) - S m ) > en q } U { inf (^(t) - S m ) < -e'n q } 

0<t<n 0<t<n 

U { sup X$ (f)+1 > (e - e')n q } 

0<t<n 

C { max \S k - ku\ > en q \ U { max \S k - ku\ > e'n q \ U { max X k > (e - 

l<fc<n l<fc<n l<fc<n+l 

where 1 < e' < e. Thus, 

exp{an qr }P( sup \n${t) - t\ > en 9 ) < exp{an 9r }P( max - fe//| > en q ) 

0<t<n l<k<n 

+ exp{an 9r }P( max \S k - kfj,\ > e'n q ) + (n + 1) enp{an qr }P(X > (e - e')n q ). 

l<k<n 

The first and the second term above converge to as n — > by Lemma A. 2. We prove 
convergence to for the third term using Markov's inequality, 

(n + l)e anqr P(X > (e - e')n q ) < (n + l) e a( - nqr -( £ - £ >^Ee aX -> as n -)• 00. 
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